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ABSTRACT 

Context. The structural parameters of the disk of the Large Magellanic Cloud (LMC) are estimated. 

Aims. We used the JH photometric data of red clump stars from the Magellanic Cloud Point Source Catalog (MCPSC) obtained from 
the InfraRed Survey Facility (IRSF) to estimate the structural parameters of the LMC disk, like the inclination, i and the position 
angle of the line of nodes (PA /o „), t/>. 

Methods. The observed LMC region is divided into several sub-regions and stars in each region are cross identified with the optically 
identified red clump stars to obtain the near infrared magnitudes. The peak values of H magnitude and (J-H) colour of the observed 
red clump distribution are obtained by fitting a profile to the distributions and also by taking the average value of magnitude and 
colour of the red clump stars in the bin with largest number. Then the dereddened peak H magnitude of the red clump stars in each 
sub-region is obtained from the peak values of H magnitude and (J-H) colour of the observed red clump distribution. The RA, Dec 
and relative distance from the center of each sub-region are converted into x, y & z Cartesian coordinates. A weighted least square 
plane fitting method is applied to this x,y,z data to estimate the structural parameters of the LMC disk. 

Results. An intrinsic (J-H) colour of 0.40 ± 0.03 mag in the IRSF SIRIUS filter system is estimated for the RC stars in the LMC and 
a reddening map based on (J-H) colour of the RC stars is presented. When the peaks of the red clump distribution were identified by 
averaging, an inclination of 25".7±1°.6 and PA/„„ = 141 .5±4°.5 were obtained. We estimate a distance modulus, yu=18.47±0.1 mag 
to the LMC. Extra-planar features which are in front as well as behind the fitted plane are identified which match with the optically 
identified extra-planar features. The bar of the LMC is found to be part of the disk within 500 pc. 

Conclusions. The estimates of the structural parameters are found to be independent of the photometric bands used for the analysis. 
The radial variation of the structural parameters are also studied. We find that the inner disk, within ~ 3°.0, is less inclined and has 
larger value of PA ;o „ when compared to the outer disk. Our estimates are compared with the literature values and the possible reasons 
for the small discrepancies found are discussed. 

Key words, (galaxies:) Magellanic Clouds; galaxies: structure; stars: horizontal-branch 

1 . Introduction The sample of RC stars used in the study by K09 have con- 

tamination from stars in the other evolutionary phases like the 
AGB stars. The details of the method applied for the reddening 

. , „ . _, , , T ,_„ S . ... , .. correction is not clearly mentioned in the paper by K09. The in- 

The Large Magellanic Cloud (LMC) is a disk galaxy with , . • ,. m , e tU r,^ . „. f , , vr , n ■ ■ 

b b , , . \ ' b . . T „„ trinsic (J -H)n colour of the RC stars estimated by K09 is in the 

planar geometry and the orientation measurements of the LMC , ■ ^ clt . ™. , . . ^. ( . 

. , • , , • , Johnson Cousins Glass niter system. This has to be tranformed 

disk plane have been estimated previously by various authors 

dde V aucouleu rs & Freeman] Il972l Ivan der Marel & Cioni 

2001, Ols en & Salvkll2002l and ISubramanian & Subramaniam 



to the SIRIUS IRSF filter system for the accurate reddening esti- 
mation. The extinction maps of the LMC estimated from the RC 

, « , , , , , , , stars data in optical bands, given in the lower panels of Fig. 5 

« 2010) using optical data of different tracers. Reddening plays an ■ c - £ c [^~r : — p c u : 1 Jorii nk a- c- 

rrt ,. ^ , • , r, , bi ^ 3 r and Fig. 6 of Subramanian & Subramaniam (2010) and in Fig. 

. . . important role in the estimation of the structural parameters of a c jrf — t ,1 L nl ,|, T r~ : r . . 

* . , rr r , , , 3 of Haschkeetal (20ll) show variation in extinction across 

galaxy. As the effect of reddening is less in longer wavelengths, , „ ... { ,, . . , 

, . . . ■ r , ,.,,„, , the galaxy. Especially a large reddening in the south-western 

the LMC structure estimated using near infrared (NIR) data ,„ ■ ■ AA ■ , •■ .-^ , 

. ,., , , , , „ , , , • i T . ' — i disk is seen in the reddening maps where warps are identified, 

is likely to have reduced effect due to reddening. Koerwer ,. . ... - / . , 

/, r t^^^s , • • • ■ r^T^ — 7 — Hence reddening correction is a important factor in the estima- 

(2009) (hereafter K09) derived an inclination, ; of 23° .5 + t . . f , , ■ Y tU , , c . , 

\ n . , . , r r , ' , tion of the structural parameters and the extra-planar features of 

4 and a postion angle of line of nodes (PA, ), of 54 .6 ±e LMQ Jn ^ anal i§ done fe KQ ^ ^ . g alsQ an oyerl 

+ l .2, using the JH photometric data of red clump (RC) stars c , ■ . ■ , , . , , . , 

r b r , _ r _ ... „ . Xi . r, ■ of sub-regions in the peripheral regions of the disk which can 

from the Infrared Survey Facility Magellanic Cloud Point . . , . „ .. . , , ^ , 

^ , „„„„ »»^r,o^,N • , cause some structural information to be averaged out. The above 

Source Catalog (IRSF MCPSC). This study was unable to .• , . . .. . , . , T . . 

, . , , , r , , , • , mentioned points motivated us to re-estimate the LMC structure 

identify the warps in the south western end of the disk, which . t , XTTT , , t c t , nr ~, . , , x^ nri T . 

. -c . , r . , ,. rrr: „ , , i u ' _ i , using the same NIR data of the RC stars used by K09. In this 

is evident in the optical studies (Olsen & Salvk 2002 and , . . . , . . • -* , TT 

' paper we use the photometnc data of the RC stars in the J and H 



Subramanian l & Subramaniam | |2010D of the structure of the bands from IRSp Mcpsc J K ato et al. Il2007t) . The K s band 
LMC using the RC stars. F ' 11 " 
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Fig. 1. A sample NIR CMD of a sub-region in the LMC is Fig. 2. A sample optical CMD of a sub-region in the LMC is 
shown. shown. The box used to identify the RC stars is also shown. 



magnitude limits of the survey is too shallow to reliably detect 
the RC stars. 

The RC stars are core helium burning stars which are metal 
rich and massive counter parts of the horizontal branch stars. 
They have a constant magnitu de which make the m standard can- 
dles for distance estimation (Stan ek et al.l[l9 98). Their constant 
characteri stic colour make the m good tracers for reddening es- 
timation (Subramaniam 2005). Here we use the JH magnitudes 
of RC stars to understand the structure of the LMC disk. The 
RC stars occupy a compact region in the optical colour magni- 
tude diagram (CMD) which makes them easily identifiable. Even 
though they occupy a compact region in the infrared CMD, their 
location is blended with the location of stars in the other evolu- 
tionary phases. In this study, we identify the RC stars in the NIR 
by cross identifying with the optically identified RC stars to min- 
imize the contamination of stars in other evolutionary phases, 
especially the AGB stars. 

The structure of the paper is as follows. In section 2 the data 
and the selection of the RC sample are explained. The analysis 
is described in section 3 and the results are given in section 4. 
The variation of structural parameters as a function of radius is 
described in section 5. The effects of reddening and population 
effects of RC stars in the estimation of the structure are discussed 
in section 6. As this work is motivated by the study of K09 and 
we use the same data set and the tracer that used in the study, 
in section 7 we compare our results with the study of K09 in 
detail. The differences seen are also explained in this section. 
The conclusions of our study are given in section 8. 




(J-H) mag 

Fig. 3. A sample NIR CMD of a sub-region in the LMC with 
the location of RC stars is shown. The red points are the RC 
stars identified from the optical CMD of the sub-region. The box 
within which the RC stars in the IR CMD are distributed is also 
shown. 



2. Data 

The IRS F Magellanic Clo ud Point Source Catalog (IRSF- 
MCPSC) (lKato etal.ll2007h is an outcome of an imaging sur- 
vey of the Magellanic Clouds (MCs) in the NIR bands J (1.25 
/mi), H (1.63 /iin) and K, (2.14 fim) during the period October 
2001 to March 2006. The observations were made with the 



SIRIUS camera (Simultaneous three colour InfraRed Imager 
for Unbiased Survey) on the InfraRed Survey Facility (IRSF) 
1.4 m telescope at Sutherland, the South African Astronomical 
Observatory. The SIRIUS camera is equipped with three 1024 
x 1024 HAWA II arrays to enable simultaneous observations i n 
the three bands (Nagashi ma et al. ll999llNagavama et al. 12 003). 
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The IRSF/SIRIUS pixel scale is 0.45 arcsec/pixel, yielding a 
field of view of 7 .7 x 7.7 arcmin 2 . The photometric catalog 
dKato et al. I l2007h includes 14811185 point sources for a 40 
deg 2 area of the LMC and 2769682 sources for an 1 1 deg 2 area 
of the SMC. In our present study, the LMC catalog is used to es- 
timate the structural parameters of the LMC using the RC stars. 

We divided the IRSF MCPSC region of the LMC into 928 
regions with a bin size of approximately 10.53 x 15 arcmin 2 . 
The average photometric error (magnitude range 15 - 20 in J and 
H bands) is around 0.1 mag. Photometric data with error less 
than 0.3 mag are considered for the analysis. A sample (/ - H) 
vs H CMD is shown in Fig. 1. To isolate the approximate RC 
location in the infrared CMD, we used the optical CMD of the 
corresponding sub-region. The RC stars are easily identifiable 
in the optical CMD as a separate component. 

The Magellanic Cloud Photometric Survey (MCPS) 
(Zaritsky et al. 1997) obtained the UBVI photometry of virtu- 
ally all stars brighter than V= 21 mag in the MCs. The five year 
survey was conducted at the Las Campanas Observatory's 1 m 
Swope telescope and the images were obtained using the Great 
Circle Camera (GCC, Zaritsky et al. 1996). The thinned 2048 
x 2048 CCD has 0.7 arcsec/pixel scale. The survey scanned 64 
deg 2 of the LMC and 16 deg 2 of the SMC. The IRSF observed 
region of the LMC comes within the MCPS observed region of 
the LMC. Zaritsky et al. (2004) presented the data of the LMC 
MCPS survey. The MCPS observed region of the LMC is also 
sub-divided into 1512 regions (each with an area of 10.53 x 
15 arcmin 2 ). For each sub-region (V-I) vs I CMD is plotted 
and the RC stars are identified. A sample optical CMD of the 
LMC is shown in Fig. 2. For all the regions, the RC stars are 
found to be located well within the box shown in the CMD, 
with boundaries 0.65 - 1.35 mag in (V-I) colour and 17.5 - 
19.5 mag in I magnitude. The average photomeric error of stars 
in the RC magnitude range is 0.05 mag in the V and I bands. 
We considered only stars with error less than 0.15 mag for the 
identification of the RC stars. 

We cross-identified the optical and infrared data and ob- 
tained the infrared (J,H) magnitudes of the RC stars identified 
within the box of the optical CMD (Fig. 2). Those are shown 
as red points in Fig. 3. Most of the RC stars are well within the 
box of infrared CMD, with boundaries 0.2 - 0.8 mag in (J - H) 
colour and 15.9 - 17.8 mag in H magnitude. Stars outside the 
box are not considered for the analysis. There are a few stars 
outside the box, especially fainter than 17.8 mag. As the H band 
limiting magnitude of the survey is 17.8 mag, it is justifiable to 
exclude those stars from the analysis. Again, from Fig. 3 we can 
see that the peak of the RC stars in H band is brighter than 17.8. 
As our method needs to identify the peak of the RC distribution, 
the box is used. The RC stars within the box are optically iden- 
tified for all the 928 sub-regions. Out of 928 regions only 926 
regions have reasonable number (100-3000) of RC stars to do 
the analysis. 



3. Analysis 

3.1. Identification of the peak magnitude and colour of the 
RC distribution 

3.1.1. Method 1 

The number distribution of the RC stars in H magnitude and 
(J - H) colour are obtained with a bin size of 0.05 mag and 




Fig. 4. A typical H magnitude distribution of a sub-region in the 
LMC. Each Gaussian profiles are shown as red and blue lines. 
The combined best fit profile of double Gaussian is shown as 
black line. The reduced ;f 2 value of the fit is 1.1. 




(J-H) mag 

Fig. 5. A typical (J - H) colour distribution of a sub-region in 
the LMC. The best fit profile to the distribution is also shown. 
The reduced^ 2 value of the fit is 1.2. 



0.04 mag respectively. The obtained distributions in colour and 
magnitude are fitted with a) a Gaussian function, b) Gaussian + 
quadratic polynomial and c) Combination of two Gaussian func- 
tions. A non linear least square method is used for fitting and the 
parameters are obtained. The parameters obtained are the coef- 
ficients of each term in the function used to fit the profile, error 
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in the estimation of each parameter and reduced x 2 value. The 
errors in the estimated parameters are calculated using the co- 
variance matrix. By comparing the reduced^- 2 values of different 
profile fits we found that combination of two Gaussian profiles, 
one narrow component and one broad component, fits well for 
the H magnitude distribution. The (J - H) colour distribution 
in majority of the (~ 60 %) sub-regions fits better with a single 
Gaussian. The colour distributions in the remaining sub-regions 
are best fitted with a double Gaussian distribution or a Gaussian 
+ quadratic polynomial. By comparing the reduced x 2 values, 
the best profile fit is used for further analysis. The H magnitude 
and (/ - H) colour distributions are shown in Fig. 4 and Fig. 
5 respectively. The best fit profiles are also shown in the fig- 
ures. From Fig. 4 we can see that the H magnitude distribution is 
well represented by two Gaussian functions. Each Gaussian pro- 
file is separately shown as blue and red lines and the combined 
profile is shown as black line. The peak of the narrow compo- 
nent coincides with the the bin with maximum number of stars. 
Thus this value is taken as the RC peak magnitude for our anal- 
ysis. The peak shift between the broad and narrow component 
ranges from 0.05 mag to 0.1 mag. The broad component peak 
is fainter than the narrow component peak in some sub-regions 
and in some other sub-regions it is brighter. The broad compo- 
nent may be representing the thick disk RC population and/or 
can be due to stars in other evolutionary phases which are the 
contaminants in the sample. The NIR sample are optically se- 
lected and in optica l magnitude distribution such a doub le peak 
feature is not seen (Subramanian & Subramaniam 2010). Hence 
we cannot conclusively say anything about the broad compo- 
nent peak and dispersion. In Fig. 5, the RC colour distribution 
which is best fitted by a single Gaussian is shown. In 20% of 
the sub-regions, the RC colour distribution shows double peak 
with broad and narrow components. The peak shift between the 
broad and narrow components is not always in one direction. 
The broad component peak is bluer than the narrow component 
peak in some sub-regions and in some other sub-regions it is red- 
der. The narrow component peak is taken as the RC peak colour. 
Note that the RC colour distribution in majority of sub-regions 
has only one peak. The parameters that are needed to estimate 
the structure of the LMC are the peak H mag, peak (J - H) mag, 
the associated errors and the reduced x 1 values. Regions with 
peak errors greater than 0. 1 mag and those with reduced x 2 value 
greater than 3.0 are omitted from the analysis. Thus the regions 
used for final analysis reduced to 775 from 926. Around 85 % of 
the sub-regions are available for the final analysis. 

3.1.2. Method 2 

In the method 1, 15 % of the sub-regions of the LMC are omit- 
ted from the final analysis due to poor fit of the RC distribution. 
Many sub-regions in the central region (within in a radius of 2.5 
degrees) are part of the omitted data points. Some of the outer 
regions are also omitted. The structural parameters of the LMC 
are very much dependent on the coverage and hence these inner 
regions w hich are omitted can make a vast d ifference in the es- 
timates ( Subramanian & Subramaniam I l2010h . The location of 
the bar with respect to the disk and the structural parameters of 
the inner region (within the radius of ~ 2.5-3 degrees) are some 
of the interesting aspects related to the structure of the LMC. 
Thus to understand the structure of the LMC in more detail, it 
is important to increase the available regions in the inner region. 
Though the most appropriate method for the estimation of the 
accurate peak of the RC distribution in a sub-region is by fitting 
a profile by numerical analysis (method 1), we can also estimate 



the peak of the distribution as the average of the magnitudes of 
stars in the bin with largest number of stars. This allows us to es- 
timate the peak mag of RC stars in all the sub-regions including 
the central regions. 

The RC stars identified are binned in both H mag and (J-H) 
colour with a bin size of 0.04 and 0.03 mag respectively. The bin 
with largest number of RC stars is identified in both the magni- 
tude and colour distribution. The average H mag, (J-H) colour 
of the stars in that particular bin and the associated standard de- 
viation are estimated. In order to reduce the effect of binning in 
the estimation of average H mag and (J-H) colour, we identi- 
fied the number of RC stars in each bin on either sides of the 
bin with largest number of RC stars. If the number of RC stars 
in those bins are greater than the number, (N - y/N) where N 
is the largest number of RC stars identified, those stars are also 
included in the estimation of average values. 

3.2. Estimation of E(J-H) reddening 

To estimate the structural parameters of the LMC disk from the 
RC magnitude we need to correct for extinction. The RC peak 
(J - H) mag at each location is used to estimate the E(J-H) 
reddening and hence the extinction in H band. The reddening is 
calculated using the relation, 

E(7 - H) = (J - H)„b s - (J - H) intrinsic 

To obtain the absolute J and H band RC magnitudes and 
hence the ( J - //),„/ n „ J/c c olour of RCstars in t he LM C, we used 
the method 1 described in lGirardi & Salarisl (1200 1 h . The mean 
RC propertie s as a function of time and metallicity based on 
iGirardi et al. I (2000) isochrones are available in table format in 
http://pleiadi.pd.astro.it We used the star formation rate from 
Harris & Zaritskv (2009) and the age metallicity relation given 
bv lPaeel & Tautvaisiene I dl998l) to obtain the absolute mean ab- 
solute J and H magnitudes of RC stars in the LMC. The RC 
stars in the LMC has an age range. They are older than 1 Gyr 
and younger than 10 Gyr. We used 1-10 Gyr age range and also 
1.5-9.5 Gyr age range and obtained the mean magnitudes. The 
mean and the standard deviation of the values obtained in the 
above mentioned two age ranges are used as the final value. The 
(J— H)i„trinsic value is obtained as 0.47 ± 0.02 mag. K09 did a 
similar analysis to estimate the intrinsic colour of RC stars, as- 
suming a constant star formation rate and a slightly different age 
range. The value obtained by K09 for (J— H) intrinsic is also 0.47 + 
0.06 mag which matches with our estimate. 

But the Girardi i sochro nes estimate the absolute magnitudes 
in iBessell & Brettl (Il988l) filters. The IRSF SIRIUS observa- 
tions are done in MKO filter system. Thus it is important to do 
the necessary transformations to obtain the (J-H) i „, r ,„ Jlc value in 
MKO system. The (J-H) i „, n „ JIC value in the MKO system has 
to be used for t he reddening cor rection. We did the necessary 
transformations ( Cai'Den tentZOO 1 1) and obtained the (J-H) i „, ri > !JIC 
colour in MKO system as 0.40 + 0.03 mag. The errors in the 
transformation co-efficients are also considered while estimating 
the final error associated with the (J— H) intrinsic colour in MKO 
system. In the process, the absolute mean H and J magnitudes 
of RC stars in the LMC are also estimated. The mean values are 
My = -1.13 + 0.02 mag and M H = -1.53 + 0.02 mag. 

The reddening, 

E(J-H) = J - H observedpeak - 0.40 ± 0.03 mag 
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towards each sub-region is estimated. The interstellar ex- 
tinction towards each sub-region is estimated using the relation, 

A H = 1.65 ± 0.16 x E(J -H) (K09). 

3.3. Estimation of the relative distance of sub-regions from 
the center 

The dereddened Ho magnitude of the RC stars for each sub- 
region is estimated using Ho = H - A#. 

The difference in Ho mag between regions is assumed 
only due to their difference in the relative distance, AD in kpc, 
calculated using the distance modulus formula given below. 

AD = 50 X l§( H aregion-H a meari)IS 

where 50 is the distance to the center of the LMC in kpc. 



3.4. Plane fitting procedure 

The R.A, Dec and relative distance of a sub-region from the 
center are used to create a cartesian coo rdinate system using the 
transformation equations given below (Ivan der Marel & Cionil 
2001, Subramanian & Subramaniam 1l2010f see also Appendix A 
of lWeinberg & Nikolaev 11200 lb . 

x = -Dsin(a - ao)cos(5, 

y = Dsini5cos(5o - Dsin5ocos(o' - ao)cosS, 

z = Do - Dsin<5shi(5o - Dcos(5ocos(o' - a^costf, 

where Do is the distance to the center of the LMC and D, 
the distance to the each sub-region is given by D = Do + AD. 
The {a, 5) and (chq , do) represents the R.A and Dec of the re- 
gion and the center of the LMC respectively. In our analysis, the 
optical center of the LMC, 05 / '19'"38'.0 -69°27'5".2 (J2000) 
dde Vaucouleurs & Freeman] (1 19721)) is taken as the center of the 
LMC. According to Ivan der Mare l & Cioni I (1200 lb the adopted 
center does not affect the derived parameters for the LMC disk 
plane (inclination, PA/ „).The distance, Do to the center of the 
LMC is taken as 50 kpc. 

Once we have the x,y, and z coordinates we can apply a 
weighted least square plane fit to obtain the structural param- 
eters of the LMC disk. The equation of the plane used for the 
plane fit is given by 

z = Ax + By + C 

From the coefficients of the plane A,B & C, i and <p can be 
calculated using the formula given below. 

Inclination, i = arccos(C/ VA 2 + B 2 + 1) 

PA /o „, <p = arctan(-A/B)+sign(B)7r/2. 

We calculated the deviations of the LMC disk from the 
plane with estimated coefficients. The expected z for a plane is 
calculated from the equation of a plane. The difference in the 
expected and calculated z values is taken as the deviation of 
the LMC disk from the plane. Thus the extra-planar features of 
the LMC disk are identified and quantified. Once the deviations 




E(J-H) mag 



Fig. 6. Two dimensional plot of reddening (E(J-H)) obtained 
based on method 1 . 




E(J-H) mag 



Fig. 7. Two dimensional plot of reddening (E(J-H)) obtained 
based on method 2. 



are estimated, the regions with deviations above the error in 
z are omitted and the plane-fitting procedure is applied to the 
remaining regions to re-estimate the structural parameters of the 
LMC disk plane. 

The uncertainty in the H and (J - H) peak values, uncer- 
tainty in the intrinsic colour of the RC stars, uncertainty in the 
reddening law are all propagated properly to estimate the error 
in relative distance. The error in the z values are used in the 
weighted least square plane fitting procedure to obtain the error 
in the estimated disk parameters. Thus the error in the estimate 
of the LMC disk parameters is calculated by propogating the 
systematic errors associated with all the quantities involved in 
the estimation. 
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4. Results 

4. 1 . NIR Reddening Map towards the LMC 

One of the by products of this study is the NIR reddening map 
towards the LMC. The shift in the peak of the (J-H) colour dis- 
tribution with respect to the characteristic (J-H) colour of RC 
stars is a measure of reddening. The peak (J-H) colour of the 
RC stars in each sub-region is estimated based on method 1 and 
also based on method 2. The intrinsic value of the E(J-H) colour 
of the RC stars in the LMC is chosen as 0.40 ± 0.03 mag. Using 
this value, the E(J-H) value of each sub-region is estimated as 
described in section 3.2. A colour coded two dimensional red- 
dening map obtained based on method 1 is shown in Fig 6. The 
E(J-H) value has a range from 0.04 mag to 0.13 mag with an 
average of 0.08 ±0.03 mag. A similar plot obtained based on 
method 2 is shown in Fig 7. The E(J-H) value has a range from 
0.05 mag to 0.15 mag with an average of 0.08±0.03 mag. From 
the figures, Fig 6 and Fig 7 we can see that there are slight vari- 
ations in the reddening estimates when region to region is con- 
sidered. But the locations where large/less reddening are seen 
coincide well. Also, the range of reddening values (0.04-0.13 
mag and 0.05-0. 15 mag) and the average reddening (0.08 ± 0.03 
mag) estimated from both the methods match well. In order to 
get a quantitative estimate of the variations, the difference in the 
reddening estimates obtained for each region using two meth- 
ods is calculated. The difference has a range from a minimum 
value of 0.0002 mag to a maximum value of 0.05 mag. In most 
of the regions, the difference is < 0.009 mag , which is much less 
than the average error (0.03 mag) associated with the estimation 
of reddening in both the methods. Only for a few regions (less 
than 6%) the difference is more than 0.03. From these maps we 
can see that reddening varies across the LMC. The reddening is 
more in the south western regions of the LMC disk as well on 
the eastern side. The reddening variations ob tained here are sim- 
il ar to those obtained in the optical stud ies dHaschke et al 11201 ll 
& ISubramanian & Subramaniam 12010b . 

4.2. Structural parameters of the LMC disk 
4.2.1. Method 1 

The structural parameters of the LMC disk are estimated using 
the dereddened mean Ho magnitude of the RC stars in 775 re- 
gions. This method gives an inclination of i= 26°. ± 1°.2 & 
PA, OB = 146°.5±3°.5 for the LMC disk. The deviation of the 
LMC disk regions from the estimated plane are calculated as 
explained in section 3.3. Deviations above 1.5 kpc are consid- 
ered as significant deviations from the fitted plane. The lower 
left and upper left panels of Fig 8 show the deviation of the 
LMC regions from the plane. In the lower left panel, all the 
regions used for the analysis are plotted. The black points are 
those which are on the fitted plane, red points are disk regions 
which are behind the plane and blue points are the disk regions 
which are in front of the plane. In the upper left panel only 
the regions with deviations above 1.5 kpc are plotted and the 
size of the points are proportional to the amplitude of the devia- 
tion. Here also, red points are regions behind the fitted plane and 
blue are in front of the fitted plane. Th e magenta points are th e 
regions where warps are identified by lOlsen & Salvk I (120021) . 
The blue and red crosses are the regions where infront of the 
plane and behind the plane features ar e identified in the study of 
ISubr amanian & Subramaniam (2010) using MCPS data. From 
the plots we can see that there are some regions which are devi- 
ated from the planar structure of the LMC disk. Even though 



the coverage of MCPS is larger than that of the IRSF MCs 
survey and that there are large number of regions omitted in 
the IR analysis, we can see that most of the locations of de- 
viations identified in the present study match with the loca- 
tions of deviations identified from the p revious optical analy- 
sis dSubraman ian & Sub ramaniam I l20i0l) . In the south western 
LMC, there are regions which are closer to us and also away 
from us. 

The presence of extra-planar features would affect our esti- 
mation of the structural parameters of the disk. Out of 775 re- 
gions, 91 regions show deviations greater than 1.5 kpc and are 
considered as real deviations. We removed these regions and re- 
estimated the structural parameters of the disk using the remain- 
ing 684 regions. We obtained an inclination of i= 26°. 6 + 1°.3 & 
PA /o „ = 148°.3+3°.8. The dereddened RC magnitude is plotted 
against the axis perpendicular to the line of nodes which is the 
axis of maximum gradient and it is shown in the lower panel of 
Fig 9. The plot clearly shows the effect of inclination from NE 
to SW of the LMC disk. The black points are those on the plane 
of the disk and the red are the extra-planar regions. The slope 
and y-intercept of the line fitted to the data, excluding the extra- 
planar regions are 0.019 ± 0.001 mag/kpc and 16.94 + 0.1 mag. 
The slope is the measure of the inclination and the inclination 
estimated from the slope is 25°.4±1°.4. 

4.2.2. Method 2 

Using this method we obtained an inclination of i= 26°. 6 ± 1°.4 
& PA /0 „ = 138°.7±3".5 for the LMC disk from the analysis of 
919 sub-regions. The standard deviations estimated for the peak 
H mag and (J - H) colour are used to estimate the error in the 
estimation of distance. The regions which show deviation larger 
than 1 .5 kpc are considered as real deviation. The lower right 
and upper right panels of Fig 8 show the deviation of the LMC 
regions from the plane. The optically identified deviations are 
also plotted. The symbols are the same as shown in the lower 
left and upper left panels of Fig. 8. Most of the deviations identi- 
fied are similar to that obtained based on method 1 . As there are 
more regions in the analysis based on method 2, there are some 
more extra-planar features identified. These deviations more or 
less match with the optically identified deviations. The structures 
in the south western region are better revealed in this analysis. 
As seen earlier, in the south western region along with the re- 
gions which are closer to us ( which matches with t he regions 
where warps are identified by lOlsen & Salvk 1 12002). there are 
regions which are away from us as well. Based on the ampli- 
tude of the deviations in the south western region, the deviations 
closer to us are significant than those away from us. Out of 919 
regions, 215 regions show deviations. These regions are removed 
and the structural parameters are re-estimated. An inclination of 
i = 25".7±1°.6 and PA,„„ = 141°.5±4°.l are obtained. The incli- 
nation obtained matches well within errors with the inclination 
obtained based on method 1 . The PA/ „ is slightly different. As 
the estimated parameters very much depend on the coverage of 
the data set, the variation in these values can be attributed to 
more number of inner and outer regions included in the analysis 
based on method 2 as compared to method 1 . The dereddened H 
magnitude is plotted against the axis of maximum gradient and 
is shown in the upper panel of Fig. 9. The colour scheme is the 
same as that of the plot in the lower panel of Fig. 9. From the 
slope and intercept of the line fitted to the data points, the incli- 
nation of the LMC disk and the distance modulus to the LMC 
center are estimated. We estimated an inclination of 25°.0±1°.5 
for the LMC disk. 
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Fig. 8. The IRSF MCPSC regions which are fitted on the plane and those which are deviated are shown. The lower left and upper 
left panels are the plots which are obtained from the analysis using method 1. The lower right and upper right panels are the plots 
obtained from the analysis using method 2. In the lower panels, black dots represent regions on the fitted LMC plane, red dots 
represent regions behind the fitted plane and the blue dots represent regions which are in front of the fitted plane. The upper panels 
show only regions with deviations, greater than 1.5 kpc out of the plane and the size of the points are proportional to the amplitude 
of the deviations. The blue and red crosses are regions which are infront o f the plane and behind the plane respectively, identified in 
th e optical study u s ing M CPS data (Subramanian & Subraman iam 112010 !). Magenta dots are regions which are suggested as warps 
bv lOlsen & Salvkl ([2002). The green hexagon in both the panels represents the optical center of the LMC. 



The significant result of this analysis, based on method 2, is 
that the structure of the central region is revealed. Most of the re- 
gions in the bar are on the plane and there are only a few regions 
which show deviations from the plane. This suggests that the bar 
is co-planar with the disk. This result is very similar to the result 
obtained bv lSubramaniam & Subramanian I (12009b regarding the 



location of the bar with respect to the disk. Even though most of 
the regions in the bar are on the plane of the LMC disk, there 
are some localized regions which are brighter than the nearby 
regions. There is also one region near the center which is away 
from us. These suggest that the bar region of the LMC is struc- 
tured. 



8 



Subramanian & Subramaniam: Structure of the Large Magellanic Cloud 




Axis of maximum gradient (kpc) 

Fig. 9. Dereddened RC magnitude plotted against the axis of 
maximum gradient. The RC magnitude estimated from method 1 
and method 2 are plotted in lower and upper panels respectively. 
The red points are regions which show deviation larger sigma. 
The direction of inclination is shown as red line. 

4.3. Distance Modulus 

The y-intercepts obtained from the Fig 9 and Fig 11 are the 
mean Ho value of the RC stars in the LMC. The mean distance 
modulus to the LMC can be obtained using the formula, 

/"0 = Homean " MfffLMQ- 

In the above equation, M# is the absolute H band magni- 
tude of the RC stars in the LMC. M H is taken as -1.53 ± 0.02 
which is calculated in section 3.2. The Ho mea „ obtained from the 
lower panel of Fig 9 (method 1) is 16.94+0. 1 mag and from the 
upper panel of Fig 9 (method 2) is 16.95+0. 1 mag. Then the 
mean distance modulus to the LMC, fiQ is 18.47+0.1 mag (from 
method 1) and 18.48+0.1 mag (from method 2). These values 
match with in errors with the previous estimates of 18.54+0.06 
mag (K09), 18.55+0 01 mag & 18.5+0.01 mag (fromMCP S and 
OGLE III data sets. I Subramanian & Su bramaniam 1l20lo|). 18.5 
± 0.02 (lAlves 112004 and 18.53 + 0.07 dSalaris & Girardi ll2002l) 
towards the LMC. 

5. Dependence of structural parameters on the 
photometric band 

In the present study to estimate the structural parameters of 
the LMC disk, we assume that the variation in the extinction 
corrected H band magnitude of the RC stars between various 
sub-regions is solely due to the distance effect induced by the 
structure of disk. Then the variations should be independent of 
the chosen photometric bands. Thus, it is useful to compare the 
results obtained from the magnitude distribution of RC stars 
in H band with that of J band analysis. The analysis in J band 
described below is done similar to method 2 given in section 
3.1.2 for H band analysis. The optically identified RC stars are 



over plotted in the J vs (J - H) CMD and is shown in Fig 10. 
From the CMD we can see that most of the RC stars are well 
within the box of infrared CMD, with boundaries 0.2 - 0.8 mag 
in (J - H) colour and 16.4 - 18.4 mag in J magnitude. For all the 
926 sub-regions, the RC stars are identified with in this box. The 
RC stars identified are binned in J magnitude with a bin size of 
0.04 mag. The bin with largest number of RC stars is identified 
in J magnitude distribution. The average J mag in that particular 
bin and the associated standard deviation are estimated. In order 
to reduce the effect of binning in the estimation of average J 
mag, we identified the number of RC stars in each bin on either 
sides of the bin with largest number of RC stars. If the number 
of RC stars in those bins are greater than the number, (N - 
y/N) where N is the largest number of RC stars identified, those 
stars are also included in the estimation of average values. The 
E(J-H) value estimated for each sub-region in section 4.1 is 
used to find the extinction in J band for that region. The formula 
used is 

A y = 2.63 ± 0.23 * E(J-H) (llndebetouw et al Il2005h 

The extinction corrected J band magnitude for all the sub- 
regions are estimated. From the extinction corrected J band mag- 
nitudes the relative distance between regions and hence the z 
co-ordinates are obtained. A weighted least square plane fitting 
procedure is applied to estimate the structural parameters of the 
LMC disk. The deviations from the plane are also estimated. The 
structural parameters are re-estimated after removing the regions 
which show deviations larger than 1.5 kpc from the plane. An 
inclination of i = 24°.4±2°.5 and PA to „ = 146°.3±6°.5 are ob- 
tained. The parameters match well with in errors with the es- 
timates from the analysis of H band data. The deviations from 
the plane are shown in Fig 1 1 . The black points are those on the 
plane. The red and blue are points are regions which are behind 
and in front of the plane respectively. The size of the red and 
black points is proportional to the amplitude of deviations. The 
deviations obtained also match well with the deviations (shown 
in Fig 8) obtained from the analysis of H band data. The dis- 
tance modulus towards the LMC is also calculated using the J 
band data and is found to be 18.47 + 0.1 mag. The results for 
the analysis of J band data match well with the results obtained 
from the analysis of H band data. 

6. Radial Variation of the structural parameters 

In this section, the variation of the structural parameters of the 
LMC disk as a function of radius is studied. As the data points 
are more in the analysis based on method 2, we make use of 
that data to find the variation of the parameters as a function of 
the radius. We estimated the parameters using the data within 
different radii, starting from 1.5 degrees to 4.5 degrees with an 
increment of 0. 1 degree. The estimates obtained are plotted in 
Fig. 12 as function of radius. The lower panel shows the vari- 
ation of PA/ on and the upper panel shows the variation of incli- 
nation. The parameters are more or less similar above the radius 
of 3 degrees. From the plot we can also see that the inner LMC, 
(within the radius of 3 degrees) is disturbed. Inclination of the 
disk increases gradually from inner to outer region till a radius 
of ~ 3.25 degrees. The PA/ on varies much within a radius of 2.8 
degrees and then remains almost constant outwards. Variations 
in the inclination and position angle of line of nodes are seen 
with in ~ 3 degree radius. Variations with in the inner region 
suggest that the inner LMC, where the bar is also located, is 
structured/disturbed. This could be due to the the effect of tidal 



Subramanian & Subramaniam: Structure of the Large Magellanic Cloud 



9 



18 - 




o 1 
(J-H) mag 






X(kpc) 



Fig. 10. A sample (J -H)vsJ CMD of a sub-region in the LMC 
with the location of RC stars is shown. The red points are the 
RC stars identified from the optical CMD of the sub-region. The 
box within which the RC stars in the IR CMD are distributed is 
also shown. 

interactions and/or mergers experienced by the LMC. After 3.25 
degrees, both the parameters of t he disk remain almost constant 
i ~ 2 6°. and PA/ „ ~ 143°.0). ISubramani an & Subramaniam I 
2010) estimated the parameters of the outer disk (regions which 
have radius greater than 3°) and that of the inner disk (data 
within in the radius of 3°). They found the inner disk to have 
lower inclination and larger PA/ „ and vice versa for the outer 
disk. As they had only a few regions in the inner 3 degree ra- 
dius they could not see the continuous trend from inner to outer 
region. Here from this analysis, we confirm that the inner struc- 
ture of the LMC is different from the outer structure and the 
inner disk is less inclined with large PA/ OB and the outer disk is 
highly inclined with less PA/ „. The increase in the inclination of 
the outer disk, which makes the north eastern part of the LMC 
more closer to us, could be due to the tidal interaction with the 
our Galaxy. 

7. Effect of reddening and population effects in the 
detection of extra-planar features 

The extra-planar features identified in the disk of the LMC are 
important as they give clues to the interaction of the LMC with 
the external systems. The effects of the population differences of 
the RC stars and reddening in the detection of the extra-planar 
features have to be looked at carefully to understand whether the 
structures identified are real or not. In this section we discuss in 
detail the above mentioned effects. 

7.1. Heterogeneous population 

The RC stars in the LMC disk are a heterogeneous popu- 
lation, and therefore they would have a range in mass, age, 
and metallicity. The density of stars in various locations will 



Fig. 11. The IRSF MCPSC regions which are fitted on the plane 
and those which are deviated based on the analysis of J band 
data are shown. The black dots represent regions on the fitted 
LMC plane, red dots represent regions behind the fitted plane 
and the blue dots represent regions which are in front of the fitted 
plane. The size of the blue and red points are proportional to the 
amplitude of the deviations. 



also vary with the local star-formation rate as a function of 
time. These factors result in a range of magnitude and color 
of the net population of RC stars in any given location and 
would contribute to the observed peaks in magnitude and 
color distributions. Therefore, the deviations found in some 
regions may also be due to these p opulation differences o f 
RC stars. But previou s studie s by lOlsen & Salvk I d2002l) 
Ivan der Marel & Cioni I (120011) . ISubramaniam & Anupama 
(l2002l) . lGrocholski et all (120061) and TPiatti et al. I (120131) suggest 
that the population effects in the magnitude and colour of the RC 
stars in the central regions of the LMC are likely to be negligible. 

The contamination of the RC sample by stars in the other 
evolutionary phases, like the AGB stars, is less in our analysis 
as the stars in our sample are optically identified RC stars, as 
described in Section 2. 

7.2. Effect of Reddening 

The extra-planar features which are found to be located both 
behind the disk and in front of the disk, could be in the plane 
of the LMC disk itself if there were an over-estimate or under - 
estim ate of the reddenin g respectively. It has been demonstrated 
bv lZaritskv et al.l d 1997b that the extinction property of the LMC 
varies both spatially and as a function of stellar population. In 
our study, the dereddening of RC stars is done using the red- 
dening values estimated from the RC stars itself. In order to un- 
derstand the effect of reddening, we plotted a two dimensional 
plot of reddening as well as the deviations. The deviations and 
the reddening values obtained by applying method 1 as well as 
method 2 are given in Fig 13. The lower left and lower right 
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Fig. 12. The variation of the structural parameters of the LMC 
disk as a function of radius is shown. The upper panel shows the 
variation of inclination and the lower panel shows the variation 
of the PA /o „. 

panels show the reddening distribution obtained from method 1 
and method 2 respectively. The upper left and upper right panels 
show the distribution of the deviations obtained from method 1 
and method 2 respectively. In the lower panels, the size of the 
point is proportional to the reddening value and in the upper 
panels the size of the point is proportional to the amplitude of 
the deviation. The red points in the upper panels of the plots rep- 
resent the regions behind the plane and blue points represent the 
regions in front of the plane. The regions in the southwestern 
part of the LMC disk around our suggested warps show more 
reddening. Thus if we assume a constant reddening across the 
LMC we will not be able to identify these warps. At the same 
time, we do not see a strong correlation between reddening and 
the deviation, because both positive and negative deviations are 
observed for regions with large reddening. That is, the reddening 
could not have been both under and over-estimated. 

8. Comparison with previous estimates 

Previously many studies have been done to obtain the pla- 
nar parameters of the LMC disk using various tracers. The 
values obtained from those studies along with our estimates 
are summarised in Table 1. Tracers used in those studies are 
also mentioned in the table. From Table 1 we can see that 
the structural parameters have a range of values. The incli- 
nation value varies from 23°±0.8° to 37°.4±2°.3 and PA /o „ 
varies from 122° .5+8°. 3 to 170°+5. Our estimates are found 
to be with in this range. Subramanian & Subramaniam (2010) 
suggested that the complicated structure of the inner LMC 
causes variation in the estimated planar parameters depend- 
ing on the area covered for each study. In the present study 
also, in section 5 we showed that the inner structure of the 
LMC is disturbed with structural variations. Vario us studies 
of the LMC disk and bar regions (Fig. 6 given in iKoerwerl 



Fig. 13. The lower left and lower right panels show the distri- 
bution of reddening across the LMC, estimated using method 
1 and method 2 respectively. The size of the points is propor- 
tional to the amplitude of the reddening. The upper left and upper 
right panels show the distribution of deviations from the plane 
estimated using method 1 and method 2 respectively. The blue 
points are the regions which are infront of the plane and the red 
points are those which are behind the plane. 

(2009),Fig.2 in Subramaniam & Subramanian I (120091) and Fig. 4 
in lSubramaniaml (120031) ) have shown that it is a highly structured 
galaxy. Thus the variations in the estimate of the planar parame- 
ters are likely to be due to these structures. 

8.1. Comparison with the study of K09 

In our present study we used the same data set and tracer used by 
K09 to understand the structure of the LMC disk. In this section 
we compare the results of their study with ours in detail. K09 
obtained an inclination of 23°.5 ± 0°.4 and PA /o „ to be 154°.6 ± 
1°.2 for the LMC disk. Our estimates of the inclination is slightly 
higher than the value obtained by K09 and we estimate a lower 
value for PA/„„. The variation in the values of PA/ on and incli- 
nation may be due the differences in the methodologies adopted 
in both the analysis. There is an overlap of sub-regions in the 
edges of the data set in the analysis of K09. In our analysis there 
is no overlap of sub-regions. The contamination of the RC sam- 
ple by the stars in other evolutionary phases is less in our study 
compared to the study of K09. Also, the reddening correction 
adopted by K09 is not very clear. If K09 adopted a reddening 
correction on a star-by-star basis, the intrinsic colour chosen by 
K09 is for the RC stars. The reddening correction on a star by 
star basis is applicable only if all the stars in the sample are RC 
stars. The sample of K09 has AGB stars and also stars in other 
evolutionary phases. A method less susceptible to contamination 
is to deredden the observed peak colour corresponding to the RC 
stars in each sub-region. Again, from Fig 1 of K09 which shows 
the location of selected sample we can see that there are stars 
which have (J-H) colour less than 0.47 mag (0.47 is the intrinsic 
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Reference Inclination, i PA;™,, d> Tracer used for the estimate 



de Vaucouleurs & Freeman (1972) 
Feitzinger e t al. (19771 
Caldwell & Coulson (1986) 
Luks&Rohlfs (1992) 
Kimetal. (1998) 
van der Marel & Cioni (2001) 
Olsen & Salvk (2002) 
Nikolaev et al. (2004) 
Persson et al. (2004) 
Koerwer (2009) 

Subramanian & Subramaniam (2010) (OGLE III data) 
Subramanian & Subramaniam (2010) (MCPS data) 
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Method 1 (H band magnitudes) 26°.6 ± 1°.3 148°.3±3°.8 Red Clump stars 

Method2 (H band magnitudes) 25°.7±1°.6 141°.5±4°.5 Red Clump stars 

Method2 (J band magnitudes) 24°.4±2°.5 146°.3±6°.5 Red Clump stars 
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Fig. 14. The regions in the bar of the LMC. The size of the red 
and blue points is proportional to the amplitude of the deviation. 



colour chosen by K09). Though we expect a spread in the colour 
(RC stars which are slightly bluer and slightly redder than the 
intrinsic colour) of the RC distribution due to population effects, 
internal reddening and photometric errors, the colour selection 
range used by K09 is very large and includes stars from other 
evolutionary phases. If K09 has adopted star by star reddening 
correction, then it is not clear how they did correction for those 
stars which are bluer than 0.47 mag. The above mentioned dif- 



Fig. 15. The regions used in the analysis based on method 1 are 
shown. The black points are regions on te fitted plane. The blue 
and red points are the regions which are infront of the plane and 
behind the plane respectively. The green points are the locations 
super giant shells given by Book et al. (2008). The cyan points 
are the location of major star forming regions 30 Dor and Nil. 
The magenta lines are the approximate location of the bar. The 
magenta hexagon is the optical center of the LMC. 



ferences could have resulted in the variation of the estimated 
parameters in the two studies. 
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Fig. 16. The regions used in the analysis based on method 2 are 
shown. The colour coding and the symbols are the same as in 
Fig 15. 



K09 found that the bar floats above the disk by 1 kpc and is 
closer to us. The bar regions are included in our analysis based 
on method 2. In our study most of the bar regions are on the 
plane of the LMC disk. As the error associated in the estimation 
of the distance is ~ 1.4 kpc, regions which show deviations less 
than 1.5 kpc (1 sigma) are not considered as the real deviations 
in our analysis. In that case if the bar is brighter than the disk 
by only 1 kpc then it is not identified as a feature closer to us. 
In order to clarify this, we plotted the regions in the bar which 
show deviation greater than 0.5 kpc along with those which are 
on the plane. The plot is shown in the lower left panel of Fig. 
14. The upper left and upper right panels show regions which 
deviate greater than 1 kpc and 1 .5 kpc respectively along with 
the regions on the plane. In all the panels, the blue points are 
those which are in front of the plane, the red are the regions 
which are behind the plane and the black points are those which 
are on the plane. From the lower panel of the plot we can see 
that in the bar all three types of regions are present. Among the 
deviated regions there are regions which are nearer to us as well 
as away from us. These suggest that the bar is part of the disk and 
is structured. The middle and the upper panel plots also suggest 
the same. As mentioned in the earlier section, the differences in 
the methodologies adopted in our study and that of the K09 may 
be the reasons for the difference in the results. 

K09 found that the major star forming region like 30 Dor 
(~ at X=-1.6, Y=0.25) and Nil (~ at X=2.27, Y=2.52) are be- 
hind and in front of the plane respectively. The super giant shell 
(given in Book et al. 2008), LMC4 (~ at X=-l.l, Y=2.45) is also 
found to be in front of the plane. In order to compare the loca- 
tions of these deviations, we over plotted the deviations iden- 
tified by K09 on the deviations identified in the present study 
based on method 1 and method 2 in Fig 15 and Fig 16 respec- 
tively. The location of 30 Dor and Nl 1 are shown as cyan points. 
In the figure we have also shown the location of super giant 
shells given by Book et al (2008) as green points. In the upper 



panel, the 30 Dor region lies exactly on top of behind the plane 
features. In the lower panel, the 30 Dor region is closer to the be- 
hind the plane feature. Thus we also find that the 30 Dor region 
is probably behind the disk. The star forming region, Nil also 
lies closer to the features located infront of the disk. The super 
giant shell LMC4 (~ at X=-l.l, Y=2.45) also lies nearer to us. 
These results match in both the studies. 

The results of K09 did not reflec t the warps iden t ified in the 
south - western regions identified bvlOlsen &Salvkl (f2002) and 
also by ISubramanian & Subramaniam I (1201 Oh from the optical 
studies. Our present study identifies the warps in the south west- 
ern region. K09 finds in this region a deviation from the plane in 
the opposite direction which is not very significant. In our anal- 
ysis also, we find some small amplitude deviations away from 
the plane in the south-western region. But along with these de- 
viations we find some regions which are significantly deviated 
from the plane in such a way that they are closer to us. The non 
identification of warps in the study of K09 may be due to two 
factors. One may be the lack of appropriate reddening correction 
by K09 for the RC stars in the south western region. From Fig. 
6 and Fig. 7 we can see that the reddening is relatively high in 
the south western regions. A large reddening was also found in 
the south western region bylQLsen & Salyk (2002). The second 
factor which may have contributed is the overlapping of sub- 
regions. In the analysis of K09 there is an overlapping of sub- 
regions in the south western end. From Fig 8 (upper panels) we 
can see that there are both in front of the plane feature as well as 
the behind the plane feature which may get averaged out if the 
regions are overlapped. 

9. Conclusions 

• The structure of the LMC disk is studied using the J and H 
magnitudes (near IR magnitudes) of the RC stars. We find an 
inclination of 26°.6 ± 1°.3 & PA/ on = 148°.3±3°.8 for the disk 
when the peaks were identified by profile fitting. When the 
peaks of the red clump distribution were identified by averaging, 
an inclination of 25°.7±1°.6 and PA /o „ = 141°.5±4".5 were 
obtained. 

• A reddening map based on (J - H) colour of RC stars is 
presented. 

• We estimate a distance modulus, /i=18.48±0.1 mag to the 
LMC. 

• The bar is found to be part of the disk, within 500 pc. 

•The inner (r< 3°.0) and outer (r> 3°) disk structures of the 
LMC are found to be different. 

• We identified extra-planar features in the disk similar to 
those identified from the optical analysis of the RC stars. 

• From the comparison of our results with the study of 
K09, who used the same data set and tracer, we found that the 
selection of the RC sample from optical identification helps 
to reduce the contamination of the sample by stars in other 
evolutionary phases. 
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